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National Aeronautics and Space Administration 
Cleveland, Ohio 

SUMMARY 

* 

The gas core reactor is a proposed device which features high tempera- 
tures. As such, it has applications in high specific impulse space missions, 
and possibly in low thermal pollution MHD power plants. In this concept the 
nuclear' fuel is a ball of uranium plasma radiating thermal photons as opposed 
w to gamma rays. This thermal energy is picked up before it reaches the solid 
cavity liner by an inflowing seeded propellant stream and convected out through 
a rocket nozzle. A ”wall-bum-out” condition will exist if there is not enough 
flow of propellant to convect the energy back into the cavity. Any real reactor 
must therefore operate with a certain amount of excess propellant flow. Due 
to the thermal inertia of the flowing propellant, the reactor can undergo power 
transients in excess of the steady- state wall burn out power for short periods 
of time. The objective of this study was to determine how long the wall burn 
out power could be exceeded without ’’burning out” the cavity liner. 

The model used in the heat- transfer calculation was one-dimensional, 
and thermal radiation was assumed to be a diffusion process. The energy 
equation was linearized about the steady state solution, and was solved (with 
a step increase in power) by using Laplace transformations. Two reactor 
. situations were analyzed; a high power (7400 MW) rocket propulsion reactor, 
and a low power (3.6 MW) cavity test reactor. 

. This study showed that the wall bum out time varies as the inverse of 
the percent increase in reactor power above cavity liner steady state bum out 
power. When wall burnout power is exceeded by only a few precent, the time 
to wall burnout is of the order of seconds for the high power reactor and 
10's of seconds for the low power reactor. If the wall burnout power is ex- 
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ceeded by about 10 percent, the wall burnout time is of the order of 0. 1 sec- 
onds for the high power reactor and seconds for the low power reactor. 

INTRODUCTION 

The gas- core reactor concept is shown in figure 1. In this concept, the 
nuclear fuel is a ball of uranium plasma radiating thermal photons as opposed 
to gamma rays. This thermal energy is picked up before it reaches the solid 
cavity liner by an inflowing hydrogen stream (propellant) which is seeded 
with sub-micron size, depleted uranium particles. (The sub- micron particles 
are used to increase the opacity of the hydrogen at low temperatures. ) The 
propellant first transpirationally cools the cavity liner and then convects the 
thermal energy out the nozzle of the reactor, (Ref. 1). This type of reactor 
has application in high specific impulse (5000 s) and moderate thrust 
(50,0.00 lb) space missions, (Ref . 2); and in low thermal pollution MHD power- 
plants, (Ref. 3). Recently a small scale, '’test reactor” modification of the 
gas-core reactor has became of interest, (Ref. 4). It consists of a gas-core 
cavity imbedded within a solid core ”driver” reactor. The purpose of this 
reactor would be to provide a relatively small size cavity to test gas- core 
reactor concepts. 

In steady state the gas core reactor heat transfer is a balance between 
radial outward diffusion of thermal photons (heat) and the radial inward con- 
vection of the propellant (coolant),, A ”wall-bumout” condition exists if 
there is not enough flow of propellant to stop the thermal photons before they 
reach the cavity liner, Ref. 1. Any real reactor will have some degree of 
fluctuation about a nominal operating power, and must therefore operate be- 
low the wall-burnout power. The farther the reactor must operate below the 
wall burnout level, the less useful the reactor; therefore, the choice of a 
nominal reactor operating power becomes a compromise between safety and 
performance. If the reactor is operating below wall-bumout power level, 
there can be power transients in excess of the wall burnout power, at least 
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for short periods of time. This is due to the ’'thermal inertia” of the hydro- 
gen propellant. The objective of this study was to determine how long the 
wall burnout power could be exceeded without "burning out”- the cavity liner. 
This time is useful as a guide for estimating the required reaction times of 
control systems for gas core reactors. 

The model used in the heat-transfer calculation was one- dimensional, 
and thermal radiation transport was assumed to be a diffusion process. The 
energy equation was linearised about the steady state solution and was solved 
(with a step increase in power) by using Laplace transformations. The 
properties of the gases and seed materials used were taken from Ref. 5-10. 

SYMBOLS 

arbitrary integration constant 
Rosseland Mean Absorption Coefficient, m“* 
arbitrary integration constant 
specific heat at constant pressure, J/ (kg-k) 
defined function of S 
inverse of F 
enthalpy, J/kg 
thermal conductivity, w/m 
Laplace transform operator 
sum index 

o 

heat flux, w/m 
Laplace transform variable 
temperature, K 
time, sec 
unit step function 
velocity, m/s 
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x length from fuel toward wall, m 

Y Laplace transform of change in enthalpy 

y integral transformation of x 

Z exponential transformation of Z 

A change from steady state operator 

3 

p density, kg/m 

a Stefan- Boltzman constant, 5. 6697x10” 8 W/(M^-K 4 ) 

6 nondimensional time 

Subscripts 

mol molecular 

o evaluated at steady state 

00 evaluated at large distance from origin 

ANALYSIS 

The model used for the heat-transfer calculation is shown in Fig. 2. 

The propellant flow enters at a relatively large distance from the origin with 
a temperature, T^. It then flows towards the origin. The positive heat flux 
from the uranium plasma is applied at the origin. The time dependent energy 
equation is written as: 


P — -pv — = - — q 
3t 3x 3x 

For a diffusion process, the heat flux is given by: 

q= _ K 9T = _JC_3h 

3x C 3x 
P 
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where 


k = M2t! + k 


3ft. 


mol 


R 


All the properties are functions of enthalpy alone, and if the enthalpy 
changes with time: 

h = h + Ah 
o 


P w p 0 + p o Ah 


K « K + K’ Ah 
o o 


The heat flux becomes: 


where 


and 


q = q Q + Aq 


= - 


K 3h 
o o 


0X 
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K 


Aq 


o 0 Ah 


C 3x 
F o 


The energy equation becomes: 


p S_Ah _ (p v ) 3_Ah _ + 0 /a o B Ah) 


0t 


0x 0x \C 0x 
p o 



6 


If we define: 



then 


3 jx a 
3y dy dx 



The energy equation can then be written as: 



8 Ah 

at 


i- f JL Ah + Ah 
9y L3y 


If the variables are transformed to 



then the energy equation becomes 


— Ah = — Ah 
96 3Z 2 


The initial condition is 


Ah(0, Z) = 0 
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The boundary condition at some large distance from the origin is 

Ah(0 s 0) = 0 

The boundary condition at the origin is 

— Ah(e,i) = /^| u(e) 

0Z \pv/ I 

(0,1) 


If 6 is Laplace transformed into S and 

Y = £(Ah) 


then using the initial condition the energy equation becomes 

a 2 

SY(S, Z) = — — Y(S, Z) 

BZ 2 

This equation can be integrated and becomes 

Y(S , Z ) = A{S)e + ^ Z + B(S)e _A ^ Z 
Applying the boundary condition at large distances from the origin 

Y(S, Q) = A + B = 0 


Applying the boundary condition at" the origin 


-A y(S,l) =(^)\ 
dZ \pv/l 

( 0 , 1 ) 


1 

s 
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and solving for A and B; the solution becomes: 

i + 's/S Z Z 


Y(S,Z)=^j 


(e H 


) 


(0,1) 


S 3 / 2 (e + ^ + e~^ 




The denominator can be expanded in a Taylor series: 


Y(S,Z) = ^j 


j jtl s -3/2f e +(Z-l-2n)'Vs _ e -(Z+l+2n)Vs'j 


( 0 , 1 ) 


If the following function is defined as: 

F(S,K) = S" 3/2 e“ K ^ 


K > 0 


then its inverse is: 


f ( e, K) 


^/ie^/^-Kerfc-^ 
Vtt 2je 


where erfc is the complementary error function. The change in enthalpy 
in response to a step increase in heat flux at the origin at time zero then 
becomes 


Ah(6,Z) 



( 0 , 1 ) 


(-i) n {f(e, 


2n + 1 - Z) - f(0,2n + 1 + Z 


>} 


The change in enthalpy in response to a pulse increase of heat flux at the 
origin and applied from time zero to time one then becomes 
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oo 

Ah(e > z)=('^9\ 'Z, (-D n 

\pv/ n=0 


f(6, 2n + 1 - Z) - f(e,2n + 1 + Z) 


( 0 , 1 ) 


- u(e - ej)[t(e - e v 2n + 1 - z) - t(e - e v 2n + 1 + z)]] 


The change in enthalpy in response to a delta function increase of heat flux 
at the origin is 


Ah(e , Z) = 



(0,1) n=0 


(-l) n £-(2n+l-Z) 2 /46> _ e -(2n+l+Z) 2 /4el 

7 t6 


The change in enthalpy in response to a arbitrary function increase of heat 
flux at the origin is 


Ah(0,Z) = 


/ e °° 

ai z 

(6\ 1) n=0 


'(-if 


7 r(e - e ? ) 


■(2n+l-zr/4(e-e ! ) 


(2n+l+Zr/4 ( 6 - 6 ')\ae' 


If it is assumed that the initial condition is steady state; then 


■(py) f 


3h a /K 9h , 

o _ _d_f 2_ 2 

3x 0x\C 0X 

\ p o 


Integrating from x to 00 gives 




K 0h 
o o 


C 3x 
F o 


= 
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Integrating from 0 to x gives 



K dh 

C p (h - h J 


where 


Q o (0) 

h JO) = -2— + h c 
(pv) r 


'O 


The calculational method used was as follows. First, the initial steady 
state temperature profile was calculated. Next, the perturbation about this 
temperature profile was calculated at the next increment in time. Third, 
the-new properties at this ’’perturbed” temperature were calculated. This 
process was repeated until the final time was reached, along with the corre- 
sponding temperature profile and associated properties. 


DISCUSSION OF RESULTS 

A cross section of the gas- core reactor is shown on Fig. 1. The pro- 
pellant flows through the cavity wall, transpirationally cooling the wall. 

Then the propellant flows from the wall toward the fissioning uranium plasma 
convecting the radiant heat that it absorbs back toward the fuel and finally 
out of the reactor. This process was modeled with a one- dimensional, time 
dependent perturbation method. The reactor described in Ref. 1; with a 
maximum power of 7400 MW and a modification of the reactor described in 
Ref. 4 with a maximum power of 3. 6 MW was analyzed. The objective of 
this study was to determine how long the maximum reactor power could be 
exceeded without burning out the cavity liner. 
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Figure 3 illustrates some of the complexity of the problem. It shows 
the "Rosseland Mean Absorption Coefficient’ 9 versus temperature. At 
relatively low temperatures, the solid seed is the dominant absorber. As 
the temperature increases, the solid- seed absorption coefficient decreases 
because the propellant density decreases - which in turn causes the number 
of solid particles per unit volume to decrease. At some temperature 
(~500Q K) the solid seed begins to vaporise, and this causes the absorption 
coefficient to decrease rather sharply, at least according to current esti- 
mates. In the temperature region between 7000-10 000 K, the absorption is 
due mainly to seed vapor, which in this case is represented by depleted 
uranium vapor properties. At temperatures above 10 000 K, the absorption 
is due mainly to the hydrogen propellant itself. 

Figure 4 shows results of some reactor power transient calculations. 
Typical radial temperature profiles in the region between the fuel and the 
cavity liner at specified times after initiation of a step power increase are 
plotted for the reactors of Refs. 1 and 4. Figure 4- A shows the reactor 
of Ref. 1. In this reactor the distance between the fuel and the cavity liner 
is about 0. 4 m. About 4. 5 kg/s of propellant flows radially inward from the 
cavity liner. The wall bum out power for this reactor is 7400 MW, and the 
initial reactor power is 9. 1 percent below the steady state wall bum out power. 
Figure 4-B shows a reactor similar to the one in Ref. 4„ In this reactor, 
the distance between the fuel and the cavity liner is about 0. 11 m and about 
0. 1 kg/s of propellant flows in through the cavity liner. The wall bum out 
power for this reactor is 3.6 MW; the initial reactor power is 9. 1 percent 
below wall bum out power. 

In Fig. 4-A the reactor of Ref. 1 is assumed to be operating at a steady 
state power that is 9. 1 percent below the cavity liner burn-out power. At 
time zero, the reactor power is given a step increase of 15 percent and is 
now operating at 105 percent of wall bum out power. Typical radial tem- 
perature profiles in the region between the fuel and the cavity liner are shown 
at times of 0, 0. 1, 0.2, 0.3, and 00 sec. At about 0. 4 seconds into the tran- 
sient period, the cavity liner bum out condition is reached. The transient 
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temperature profiles are similar to the steady state profiles. A tempera- 
ture ’’front" seems to propagate towards the cavity wall with diminishing 
velocity during the power transient,, 

In Fig. 4-B the reactor similar to Ref. 4 is assumed to be operating at 
a steady state power that is 9. 1 percent below the cavity liner burn-out power. 
At time zero, the reactor power is given a step increase of 15 percent and 
is now operating at 105 percent of wall burn out power. Typical radial tem- 
perature profiles in the region between the fuel and the cavity liner are 
shown at times of 0, 1, 2, 3, 4, 5, and °° sec. At about 5 seconds into the 
transient period, the cavity liner burn-out condition is reached. The tran- 
sient temperature profiles are again similar to the steady state profiles. 

The bum out time is of course a function of the degree to which the 
transient power exceeds the wall bum-out power. Figure 5 shows the time 
to wall bum out versus the percent power increase above initial steady state. 
The initial steady state power is either 9.1 or 16. 7 percent below wall bum 
out. That is an increase of initial steady state power of either 10 or 20 per- 
cent will result in wall bum out at infinite time. Figure 5- A shows the cases 
for the reactor of Ref. 1, and Fig. 5-B shows the case for the reactor simi- 
lar to Ref. 4. 

The time to wall bum out is longer for less severe power transients. 

For example a 25 percent increase in initial steady state power results in 
wall burn out at 0. 15 sec for the high power reactor operating at 9. 1 percent 
of wall burnout power; while for the same reactor and transient operating at 
16. 7 percent of wall bum out power results in wall bum out at something 
greater than 1 sec. The same thing is seen for the low power reactor except 
the time to wall bum out is about ten times as long. The conclusion that can 
be reached is that when wall burnout power is exceeded by only a few percent 
the time to wall burn out can be of the order of seconds Tor the high power 
reactor and 10 of seconds for the low power reactor. If the wall burnout 
power is exceeded by greater then 10 percent, the wall bum out time is the 
order of 0. 1 second for the high power reactor and seconds for the low power 
reactor. The low power reactor can probably be contiuously controlled with 
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power pulse in excess of wall burnout, while the high power reactor can not. 
The high power reactor will have to be designed not to exceed 110 percent of 
wall burnout. If it does a scram will have to be initialed and the reactor will 
have to shut down in tenths of seconds. 

CONCLUSIONS 

The temperature response of a gas- core reactor was analysed when sub- 
jected to a power transient. The objective of this study was to determine 
how long the maximum steady state reactor power could be exceeded without 
burning out the cavity liner. The relatively large magnitudes of reactor 
power transients considered here would probably be precursors of a reactor 
scram, rather than normal reactor fluctuations. 

The wall burn-out time varies as the inverse of the percent increase in 
reactor power above cavity liner steady state burnout power. The transient 
temperature profiles appear to be similar to steady state temperature pro- 
files. When wall burnout power is exceeded by only a few percent, the time 
to wall burnout can be of the order of seconds for the high power reactor 
and 10 ? s of seconds for the low power reactor. If the wall burnout power is 
exceeded by about 10 percent, the wall burnout time is of the order of 

0. 1 . seconds for the high power reactor and seconds for the low power reactor. 
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tb) LOW POWER REACTOR, STEADY STATE WALL BURN OUT 
POWER IS 5.6 MW. (INITIAL STEADY STATE POWER 
EQUALS 3. 27 MW. ) 


Figure 4. - Typical radial temperature profile in the region 
between the fuel and the cavity liner at specified times 
after initiation of a 15 percent step power increase. 




(a) HIGH POWER REACTOR, (STEADY STATE 
WALL BURN OUT POWER EQUALS 7400 MW). 



STEADY STATE POWER 

(b) LOW POWER REACTOR (STEADY STATE WALL 
BURN OUT POWER EONAIS ^ r y,V7). 

Eicjoii. D. - Time to cavity wall burnout during power 
surge. 
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